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Development and Testing of 
an Axial Halbach Magnetic Bearing 

Dennis J. Eichenberg, Christopher A. Gallo, and William K. Thompson 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

The NASA John H. Glenn Research Center has developed and tested a revolutionary Axial Halbach 
Magnetic Bearing. The objective of this work is to develop a viable non-contact magnetic thrust bearing 
utilizing Halbach arrays for all-electric flight, and many other applications. This concept will help to 
reduce harmful emissions, reduce the Nation’s dependence on fossil fuels and mitigate many of the 
concerns and limitations encountered in conventional axial bearings such as bearing wear, leaks, seals and 
friction loss. The Axial Halbach Magnetic Bearing is inherently stable and requires no active feedback 
control system or superconductivity as required in many magnetic bearing designs. The Axial Halbach 
Magnetic Bearing is useful for very high speed applications where conventional bearings cannot be used 
including turbines, instrumentation, and medical applications. In addition, this technology has potential 
application in ultra-efficient motors, computer memory systems, manufacturing equipment and space 
power systems such as flywheels. 


Introduction 

The NASA Glenn Research Center has a wealth of experience in Halbach bearing technology through 
the Fundamental Aeronautics Program. The goals of the program include improving aircraft efficiency, 
reliability, and safety. The Axial Halbach Magnetic Bearing was developed under the Fundamental 
Aeronautics Program. 

The electromagnetic concept uses permanent magnet elements attached to the side of the rotor, and 
wire coils placed in the stator shell. The permanent magnets are arranged in a “Halbach” configuration 
which results in the production of a sinusoidally varying, periodic magnetic field in the vicinity of the 
stator coils. This magnetic array configuration was pioneered by Klaus Halbach for use in particle 
accelerators. When set in motion, the time varying magnetic fields interact with the passive coils in the 
stator assembly to produce repulsive forces between the stator and the rotor providing magnetic 
suspension. The advantage of this technique is that it is inherently stable once the rotor reaches a critical 
speed, and thus requires no active feedback control or superconductivity as seen in many traditional 
implementations of magnetic suspension. 

Theoretical derivations have been developed to predict the levitation forces generated by a circular 
Halbach array and coil assembly. Finite element analyses were performed to validate the theoretical 
derivations. Experimental hardware was successfully designed and developed which served to validate 
the basic principles described and the theoretical work that was performed. 

In addition to aircraft engines, this technology has potential application in ultra-efficient motors, 
computer memory systems, instrumentation systems, medical systems, manufacturing equipment, and 
space power systems, such as generators and flywheels. 
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Electromagnetic Analysis of The Axial Halbach Magnetic Bearing 

Magnetostatic Field 

Figure 1 shows a typical arrangement of magnets comprising an axial Halbach array. In this case 
there are N m = 32 sector-shaped permanent magnets with inner radius r u outer radius r 2 and axial 
thickness T. There are four magnets per Halbach wavelength in the azimuthal (i.e.,$) direction. Each 
magnet has an index, 5 = [0, 1, ... 7V m -l]. The magnets each have a magnetization M = ±4nB r /ju 0 whose 
direction is indicated by the arrows. B r is the remanent magnetization of the permanent magnet material. 
Consider two overlapping coordinate systems, one Cartesian and the other cylindrical. The cylindrical 
plane aligns with the Cartesian x-y plane, so that the z coordinates of each system are identical. The +x 
axis of the Cartesian system lies at $ = 0 in the cylindrical system. In both coordinate systems the +z axis 
corresponds to the axis of rotation and the flat bottom faces of the sectors lie at axial position z\ = 0 with 
z 2 = T. The 5 = 0 magnet is magnetized in the axial (+z) direction. The +x axis bisects the flat bottom 
surface of this magnet. Magnets lying on the -x, +y, and -y axes are also magnetized in this direction. The 
5th magnet has a central radial axis lying at angle jS s = In s/N m . Each magnet subtends an angle in the (j) 
direction of - P\ = 2n/N m radians, where and J3\ are the locations of the side faces. 

The magnetic field component solutions at an arbitrary point in space are expressed in cylindrical 
coordinates as 


B{r,j,z) = {B r ,B t ,B,} 


( 1 ) 


with an equivalent Cartesian representation of 


B{x , y, z ) = {. B r cos (j) - B f) sin B r sin $ + B (j) cos B z } 


( 2 ) 



Figure 1. — The axial Halbach array for N m = 32 magnets. Arrows indicate the direction of magnetization 
for each individual magnet. This particular arrangement will concentrate the B field below the ring and 
cancel it above the ring. 


NAS A/TM— 2006-2 14357 


2 


Using the principle of linear superposition, the aggregate field components become the sum of 
individual contributions from equivalent surface currents on the faces of each sector in the array. This 
method of analysis derives the B field as the curl of the vector potential, as performed in (ref. 1). The 

expressions for the field components in the cylindrical system at a spatial location r = {rr + (/)(/) + zz) 
become: 
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( 5 ) 


These functions do not necessarily have closed form solutions, but their numerical implementation is 
relatively straightforward in a variety of commercial mathematical analysis software packages. 
Mathematica v5.2 provided the results for this paper. 

As demonstrated elsewhere (ref. 1), at sufficiently large values of gap distance g, each of the three 
field components exhibits sinusoidal behavior versus azimuthal position, $ Additional harmonics in the 
field components arise at axial distances that are small compared to the magnet thickness. However, for a 
large enough stator winding these harmonics are smoothed out in the flux calculations. All three field 
components also exhibit decaying exponential behavior with increasing gap distance, g. 

It would be desirable to arrive at a closed form approximation for the dynamic behavior of the 
rotating array. However, as we shall see, the dependence of the field components as well as the Halbach 
wave number and the stator winding factor on radial position complicates this approach. It is therefore 
necessary to compute the flux numerically through a winding of known spatial position and properties. 

Note that both the B^r) and B z (r) field components are available for linking flux in the stator 
winding. An r-z flux plane links the B $ field, and an r-(j) flux plane links the B z field. The following 
analyses consider both choices for the orientation of the stator winding. 

Flux, EMF, Current and Force Calculations for an r-z Stator Pole Piece 

First, position the axial Halbach array just above a stator pole piece which is sector shaped and wound 
vertically, i.e., the plane of the center of the stator winding lies in the r-z plane at azimuthal location^ = 0. 
The normal vector of this plane coincides with the (/> direction. The winding consists of N t turns that are 
wound as radial spokes pointing toward the central axis of rotation of the array. The width of the winding 
in the radial direction is W and the height in the axial direction is H. W is chosen to be approximately the 
magnet width. H is chosen so as to optimally trade off between making the Lorentz force interaction 
between the returning current and the B $ field as small as practical while minimizing the stator resistance, 
which tends to reduce the lifting force and increase the drag force and power dissipation. A gap, g, exists 
between the bottom surface of the magnets and the center of the top conductors of the winding. The axial 
location of the top of the winding is Z 2 = -g and the bottom is located at Z\ = -(g + H). Figure 2 
illustrates these concepts from a view located below the winding. 
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Figure 2. — The r-z stator winding viewed from below 
the magnet array. The azimuthal span of the winding 
matches the width of one magnet. H and W refer to 
the winding height and width, respectively. A gap, g, 
(not shown) exists between the conductors of the 
winding and the magnet surface. 


The stator flux is calculated numerically by integrating the field component normal to the winding 
plane over the surface enclosed by the stator winding. For an r-z winding the flux as a function of angular 
position is given by 


Z 2 R 2 

O r _ z (c/))=N t w^ B ( ^(r)^dS = N t w ^ ^ B ( / ) (r,<f),z)drdz (6) 

S Zi R x 

using Z 2 = -g and Z\ = -(g + H) and W = R 2 - R\ = r 2 -r\. The winding factor, w =wiw 2 , accounts for the 
spatial variance of the B field over the azimuthal span and the axial depth of the winding, respectively. 
Two properties of the B field determine the two winding factor components. First, only the fundamental 
sinusoidal component of the field versus azimuthal position will contribute to force production. The first 
winding factor component w\ then becomes the average over all winding turns when the peak of a cosine 
function coincides with the center of the winding. Secondly, the B field decays exponentially with axial 
position according to the Halbach wave number, k (refs. 1 and 2). The second winding factor, w 2 , then 
becomes the average over all layers of this axial decay factor when the top layer is normalized to unity. 
For a stator winding of N tz layers with N t< /, turns in each layer, there are N t = N tz N t( f, total turns. The first 
winding factor component w\ is therefore given by: 
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and the second component w 2 is given by 
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w 2 (r) = Te 
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( 8 ) 


tz 72=0 


where d is the diameter of the conductor comprising the stator winding, s is the average distance between 
the individual turns, and R x is the inner radial position of the winding. Unlike a linear Halbach array 
(ref. 2) or even a radial Halbach array, the axial Halbach array presents a unique situation because k is not 
a constant but a function of radial position given by, 

k(r)=N m /4r (9) 


If the winding spans approximately the same angle as one magnet and N = l layer, then w - 0.9 and 
this factor may be taken outside the integral. For other winding geometries the winding factor may be 
computed as a spatial average over the radial span and axial depth of the winding and then taken outside 
the integral. 

Using the sinusoidal nature of the flux versus azimuthal position to solve equation (1) numerically for 
the maximum flux (at an angular position 2n/Nf) gives an approximation to the flux as 
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( 10 ) 


Assume that the array now rotates at a constant angular velocity G)=27rx rpm/60, where rpm is the 
rotational speed expressed in revolutions per minute. The initial rotary position of the array relative to the 
static Halbach array case is Therefore the spatial component (j) of the vector r (r, (j),z) is now a function 
of time such that 


(/) = (/). -cot (11) 

and the field components therefore become sinusoidal functions of time 


(r, -cot, z ) = -B^ o ( r)e~ k(r)z sin(®7) 

(12) 

B z (r, -cot, z) = B z0 (r)e~ k{r)z cos {cot) 

(13) 
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for simplicity, choose (j) t = 0, so that the central axis of the +z polarized magnet whose index is 5 = 0 
corresponds to the Cartesian x-axis (</>= 0). Since there are NJ4 Halbach arrays in the disk, the electrical 
frequency becomes N m co/4. Using equation (5) and Faraday’s law, the open circuit voltage becomes 
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whose peak value at time t = 0 is 
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(14) 


(15) 


computing the maximum current from equation (10) and the complex impedance of the winding via 
Ohm’s law gives 
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with R tot and L tot representing the total resistance and inductance, respectively, of the stator winding 
circuit, including the stator self-resistance, R s , self-inductance, L s , plus any external passive components. 
Typically, an external coil with resistance R c and inductance L c is added in series with the stator winding. 
The purpose of the coil is to optimize the phase shift of the voltage and current, which maximizes axial 
force production (lift) while minimizing transverse plane force production (drag). For the arrangement 
under consideration, it is necessary that the phase angle between the voltage and current be as close to n/2 
as practical so that the current is flowing radially outward (+r) when a polarized magnet is passing 
over the stator so that lift force (-z) may be produced. Adding the coil resistance and inductance gives R tot 
= R s + R c and L tot = L s + L c which gives the amplitude of the winding current as 
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The current may be decomposed into its imaginary (lift-producing) and real (drag-producing) 
components: 
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(19) 
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Note that the lift current is in phase with the inverse of the B $ field component given in equation (7) 
and the drag current is in phase with the B z component given in equation (8). In order to approximate the 
Lorentz force on the stator winding using dF = Idl x B, one first requires the spatial averages of the peak 
field components over the radial span of the top of the winding. 

-1 Rl 

„ f B^r,cot = nl2,z = -g)dr (20) 

z F /?, 

1 r 2 

= -T—r J BAr, cot = 0, z = -g)dr (2 1 ) 

These expressions have straightforward numerical solutions. B $ is a sinusoidal function with 
amplitude B^g^ and B z is a co sinusoidal function with amplitude B zavg>top . The time average of the 
product of two in-phase sinusoids is one-half the product of their amplitudes, so the Lorentz force 
becomes 


N, ™l u , ( 22 ) 

1 ~^H 

=^ B „ p N t wI Dmg W( l-e ^ ) (23) 


The exponential factor at the end of each expression accounts for the bottom leg of the stator winding, 
where the current flows in the opposite direction and force production opposes the contribution from the 
top leg of the winding. From the previous expressions, it should be noted that both the lift and drag forces 
increase with the square of the magnet strength (B r 2 ), the square of the number of stator turns (N t 2 ) and the 
square of the winding factor (w 2 ). Force production falls off with gap distance g by a factor of e 2g . The 
stator winding circuit resistances and inductances determine both the saturation value of the lift force and 
the rotational speed at which saturation is achieved. The drag force outpaces the lift force at low rotational 
speeds but reaches a maximum value at a critical speed 


CO L =- 


4 K+K 

s c_ 

N L +L 


(24) 


and then rolls off toward zero at high speeds. The expressions for lift and drag force are per stator pole, so 
the net forces must be summed over all poles. 

Flux, EMF, Current and Force Calculations for an Stator 

For the r-(j) stator the analysis proceeds in a similar way. Consider a stator pole piece which is sector- 
shaped and wound horizontally, i.e., the plane of the center of the stator winding lies in the r-$ plane 
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centered at azimuthal location $5=0. The normal vector of this plane coincides with the z direction. The 
winding consists of N t turns and the width of the winding in the radial direction is W The winding spans 
an angle J3 in the azimuthal direction. W is again chosen to be approximately the magnet width. /? is 
chosen to be exactly the angular span of two magnets or /? = 47r/Nrn so that each radial leg of the winding 
will carry a current that is in phase with the B $ component to produce lift. A gap, g, exists between the 
bottom surface of the magnets and the center of the top conductors of the winding. Figure 3 illustrates 
these concepts from a view located below the winding. 

The stator flux for this configuration is 

In 

Ri K, 

^ r -^</>)=N t wjB z (r)»dS = N t w J | Bf(r,0,z)rd0dr (25) 

S R x 2n 

N m 

The winding factor, w = W\W 2 , for a layer of N t = N tz N t( /, turns is still given by equations (3) and (4). 
The approximation of the stator flux becomes 



Figure 3. — The r-cf> stator winding viewed from below the magnet array. The azimuthal 
span of the winding (3 matches the width of two magnets W refers to the winding width. 
A gap, g, (not shown) exists between the conductors of the winding and the magnet 
surface. 
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and the open circuit voltage becomes 
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The voltage reaches a peak positive value at time t = 6n:/N m (o given by 
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The maximum current is found from the magnitude of the complex impedance as before 
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(29) 


where again R tot = R s + R c and L tot = L s + L c to account for any external inductance required to produce the 
desired phase shift between voltage and current. Interestingly, the r-(f) stator has an inherent -n/2 phase 
shift arising from the geometry between the center of the winding and the current-carrying outer legs. 
However, this is negated by the inherent n/2 phase shift between the field producing the flux ( B z ) and the 
field driving the lifting force ( B ^). Thus, it is still necessary that the phase angle between the voltage and 
current is as close to n/2 as practical so that the current is at its maximum value when ±0 polarized 
magnets are passing over the two radial legs of the stator winding. 

The current may again be decomposed into real (lift-producing) and imaginary (drag-producing) 
components: 
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(31) 


Although the flux is now produced by the B z component the lift force is still produced by the B $ 
component. So the lift current is still in phase with the B $ field component and the drag current is in phase 
with the B z component. Unlike the r-z stator, the r-(j) stator produces lift and drag forces in the same 
direction and of equal magnitude from both legs of the stator winding. This yields the following 
expressions for the average lift and drag forces. 
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( 32 ) 
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As with the r-z stator case, both the lift and drag forces increase with the square of the magnet 
strength ( B r 2 ), the square of the number of stator turns (N, 2 ) and the square of the winding factor (w 2 ). 
Force production falls off with gap distance g by a factor of e~ 2§ . The stator winding circuit resistances 
and inductances determine both the saturation value of the lift force and the rotational speed at which 
saturation is achieved. The drag force outpaces the lift force at low rotational speeds but reaches a 
maximum value at a critical speed 


(D„ 


4 K+K 

N m L s +L c 


(34) 


and then rolls off toward zero at high speeds. 

Comparative analysis of the two stator designs reveals that the force production per pole more than 
quadruples when the r-(j) stator design is used. However, packaging constraints will permit only half as 
many r-(j) stator poles as r-z poles. The net gain in force is therefore slightly more than doubled. However, 
this applies to both lift and drag forces, as well as to power dissipation. The gain in levitating force comes 
with no significant gain in lift/drag ratio and at the cost of increased real power dissipation and heating. 
The choice of stator design represents a tradeoff among these factors. 


Test Objectives 

The objective of testing the Axial Halbach Magnetic Bearing is to validate the derived theoretical 
analyses and associated finite element analyses. This validation provides confidence in using the derived 
and numerical analyses for developing and optimizing conceptual Axial Halbach Magnetic Bearing 
designs. 

Testing of the Axial Halbach Magnetic Bearing was performed at NASA Glenn. Of particular interest 
is the performance of the bearing. Force, voltage, current, and temperature were monitored at various 
rotor operating speeds. 


Test Hardware Description 

The Axial Halbach Magnetic Bearing Test Model hardware is shown in figure 4 and a block diagram 
of the system is shown in figure 5. The Test Model components are identified in figure 6. This test 
hardware includes the support structure, force measurement scale and positioning system, and drive 
system components. 

The support structure consists of precision modular structural components that are securely fastened 
to a heavy optics plate. The motor and the drive shaft bearing blocks are connected to right angle brackets 
that are connected to the optics plate via posts. This subsystem supports and aligns the drive system with 
respect to the force scale and positioning system that is placed on the optics plate. A safety shield is also 
incorporated into the support structure to protect personnel from the rotating components of the drive 
system and rotor. 

The drive system consists of a variable speed computer-controlled servo motor that is connected to a 
one half inch diameter drive shaft through a flexible coupling to allow for misalignment. The drive shaft 
is supported by two roller bearings that are secured in the bearing blocks. The rotor assembly is pinned to 
the drive shaft and is driven by the motor to maintain precise operating speed. 
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Figure 5. — Axial Halbach magnetic bearing test model block diagram. 
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The rotor assembly, as shown in figures 7 and 8, consists of a 4.7 in. overall diameter aluminum hub 
that is bored out to contain 32 Neodymium Iron Boron (Nd) B55 magnets. The magnets are each 1 quarter 
in. thick by 1 in. long segments. The 32 segments comprise a ring with an outside diameter of 4 in. and 
inside diameter of 2 in. Each magnet segment occupies 1 1.25° of the assembled magnet ring. The 
magnets are covered by a total of eight 45° 300 series stainless steel segmented covers with each cover 
restraining four magnets. This segmented cover design was necessary due to the difficulty in handling and 
assembling the 32 magnets into a Halbach array since adjacent magnets repel each other. Four magnets 
were oriented in a Halbach array and held in position and the cover was fastened to the hub with non- 
magnetic stainless steel screws. The total mass of the assembled rotor is 660 g (1.46 lb). 
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The rotation of the rotor causes the magnetic field to react with a coil of wire wound on the stator 
pole piece. Four different stator pole pieces were tested. Two r-z pole pieces were developed, one with 
one turn of wire, and the other with six turns of wire. Two r-(/) pole pieces were also developed, one with 
one turn of wire and the other with six turns of wire. An r-z pole piece is shown in figure 9 and an r-cj) 
pole piece is shown in figure 10. All of the pole pieces have 20 gauge square copper magnet wire. The 
wire is wound around a slot in the pole piece machined from Ultem polyetherimide with a maximum 
operating temperature of 340 °F. The magnet wire is cemented in place in the slot to ensure that accurate 
data is obtained. 



Figure 9. — Axial Halbach magnetic bearing test model r-z stator 
pole piece. 



W 


Figure 10. — Axial Halbach magnetic bearing r-(|) stator pole piece. 
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The pole piece may be terminated by an inductor to introduce a phase shift to maximize the generated 
repulsive force. The inductor consists of 17 turns of 10 gauge square copper wire wound around a ferrite 
rod. 

The orientation of the pole piece is highly critical. The pole piece is connected with nylon fasteners to 
a precision horizontal linear stage which is then connected to a precision vertical stage. These two stages 
allow for an accurate adjustment of the coil placement in close proximity to the surface of the rotor. The 
rotation of the Halbach array adjacent to the pole piece assembly generates a repulsive downward force 
on the force measurement scale. 


Instrumentation 

The Axial Halbach Magnetic Bearing Test Model was instrumented to measure bearing performance. 
Data from the drive motor, particularly rotor speed, were obtained from the motor drive system and sent 
directly to a personal computer (PC), where the data was stored. The other data were sent to the various 
test instruments and monitored. Test data included stator pole piece voltage, stator pole piece current, 
stator pole piece force, stator pole piece temperature, and ambient temperature. 

A block diagram of the instrumentation system is shown in figure 1 1 . Type K thermocouples were 
used for all temperature measurements. Hall Effect transducers were used for all current measurements. 



Figure 11. — Axial Halbach magnetic bearing instrumentation block diagram. 
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Test Procedures 


The tests described in this report were conducted at the NASA Glenn Research Center in Cleveland, 
Ohio. The tests were conducted in accordance with the test matrix provided in table 1. 

TABLE L— AXIAL HALBACH MAGNETIC BEARING TEST MATRIX 


Parameter 

Gap 

(in.) 

Speed 

(rpm) 

Termination 

Six turn r-z stator pole 

0.05,0.08,0.11 

100, 200, 500, 1000, 2000, 
3000 

Open circuit, short circuit, 
inductor 

Six turn r-z stator pole 

0.05,0.08,0.11 

100, 200, 500, 1000, 2000, 
3000 

Open circuit, short circuit, 
inductor 

Single turn r-Q stator pole 

0.05,0.08,0.11 

100, 200, 500, 1000, 2000, 
3000 

Open circuit, short circuit, 
inductor 

Six turn stator pole 

0.05,0.08,0.11 

100, 200, 500, 1000, 2000, 
3000 

Open circuit, short circuit, 
inductor 


The test procedure used four different stator pole designs: single-turn r-z plane, 6-tums r-z plane, 
single turn r -(/) plane and 6-tums r-^ plane. All four stators contain 20 AWG wire ( d = 0.032 in. 

(0.8 mm)). The test consisted of precisely locating the stator beneath the magnet array such that the air 
gap was one of three selected values: 0.05 in. (1.3 mm), 0.08 in. (2.0 mm, and 0.1 1 in. (2.8 mm). 
Factoring in the thickness of the retaining cover plate and the finite radius of the conductor, the effective 
gap distances from the magnet surfaces to the center of the conductor were 0.096 in. (2.4 mm), 0.126 in. 
(3.2 mm, and 0.156 in. (4.0 mm), respectively. With the array rotating at fixed mechanical speeds of 100, 
200, 500, 1000, 2000, and 3000 rpm, peak-to-peak open circuit voltage (V oc ) and short circuit current 
(Iss) were measured via an oscilloscope. With a series coil (17 turns, 10 AWG, L c = 9.1 /z/f, R c = 2.5 m/2) 
inserted into the stator winding circuit, current (I c ) and lift force (Fi ift ) production were measured with a 
precision scale. 

The analytical model presented in this article predicted expected values for V oc , I ss , I c , and Fi ift . 
Predicted values were compared to the measured results. Additionally, the analytical model predicted the 
drag force (Fdrag) and power dissipation (Pd). 

Finite element analysis (FEA) predictions at 500 and 2000 rpm at an air gap of 0.05 in. were also 
compared to the analytical and measured results for open-circuit voltage and short-circuit current. 
Maxwell 3D (Ansoft Corp.) was the software used to generate the FEA predictions. 

Test Results 

System Performance 

Various tests were conducted to determine system performance, per table 1. Table 2 summarizes the 
results obtained at 500 rpm and table 3 summarizes the results obtained at 2000 rpm for the four different 
stator designs. A complete set of data for all rotational speeds is presented as a series of plots in appendix B. 


TABLE 2.— AXIAL HALBACH MAGNETIC BEARING PERFORMANCE AT 500 rpm 


At 500 rpm 
eff gap = 0.096 in. 
(2.5 mm) 

1 Turn r-z stator 

6 Turn r-z stator 

1 Turn r-(j) stator 

6 Turn stator 

Quantity 

Measured 

Analytical 

predicted 

FEA 

predicted 

Measured 

Analytical 

predicted 

FEA 

predicted 

Measured 

Analytical 

predicted 

FEA 

predicted 

Measured 

Analytical 

predicted 

FEA 

predicted 

Voc(V) 

0.046 

0.047 

0.046 

0.231 

0.249 

0.260 

0.095 

0.100 

0.12 

0.562 

0.588 

0.7 

Iss (A) 

6.2 

6.0 

5.9 

11.8 

11.1 

11.7 

8.8 

9.9 

11.0 

24.7 

26.2 

32.0 

foil (A) 

4.6 

4.4 

N/A 

9.8 

10.0 

N/A 

7.2 

7.6 

N/A 

27.0 

30.2 

N/A 

Fiift (N) with coil 

N/A 

0.005 

N/A 

0.028 

0.028 

N/A 

0.013 

0.014 

N/A 

0.172 

0.18 

N/A 

Fdrag (N) with coil 

N/A 

0.013 

N/A 

N/A 

0.168 

N/A 

N/A 

0.045 

N/A 

N/A 

0.91 

N/A 

Pd (W) with coil 

N/A 

0.042 

N/A 

N/A 

0.582 

N/A 

N/A 

0.167 

N/A 

N/A 

3.59 

N/A 
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TABLE 3.— AXIAL HALBACH MAGNETIC BEARING PERFORMANCE AT 2000 rpm 


At 2000 rpm 
eff gap = 0.096 in. 
(2.5 mm) 

1 Turn r-z stator 

6 Turn r-z stator 

1 Turn r-(j) stator 

6 Turn stator 

Quantity 

Measured 

Analytical 

predicted 

FEA 

predicted 

Measured 

Analytical 

predicted 

FEA 

predicted 

Measured 

Analytical 

predicted 

FEA 

predicted 

Measured 

Analytical 

predicted 

FEA 

predicted 

V 0C (V) 

0.18 

0.19 

0.175 

0.94 

1.00 

0.96 

0.36 

0.40 

0.49 

2.22 

2.35 

3.00 

Iss (A) 

23.4 

24.8 

23.3 

42.20 

43.86 

42.0 

34.00 

39.98 

49.00 

105.00 

118.74 

144.0 

I coil (A) 

10.00 

10.22 

N/A 

32.50 

32.48 

N/A 

19.00 

20.17 

N/A 

77.50 

80.79 

N/A 

Fii ft (N) with coil 

0.03 

0.03 

N/A 

0.27 

0.30 

N/A 

0.09 

0.09 

N/A 

1.24 

1.74 

N/A 

Fdrag (N) with coil 

N/A 

0.02 

N/A 

N/A 

0.46 

N/A 

N/A 

0.08 

N/A 

N/A 

2.17 

N/A 

P d (W) with coil 

N/A 

0.00 

N/A 

N/A 

4.43 

N/A 

N/A 

0.51 

N/A 

N/A 

20.26 

N/A 


Discussion 

The objective was to verify the analytical model of the axial Halbach array presented in this article by 
comparing both measured results and FEA-predicted results. The analytical model effectively predicts 
induced EMF, current and lifting force in the stator pole at lower speeds. At higher speeds, the combined 
effects of temperature rise and skin effect (both of which tend to increase the effective resistance of the 
stator winding) attenuate the current and lifting force. Saturation of the ferrite core in the external series 
inductor limits the amount of phase shift that can be achieved above around 50 amps of total current. This 
also contributes to losses in force production. 

The r-(j) stator produces about 5 times the net lifting force compared to the r-z stator. This is due to 
three effects. First, the B z field is slightly stronger (-2 percent) than the B $ field. Second, the proximity of 
the entire flux plane to the magnets in the case produces more efficient flux linkage in the stator 
winding. Finally, both legs of the stator winding contribute to force production harmoniously in the r-(j) 
case, whereas the bottom leg of the r-z winding produces forces in opposition to the forces produced by 
the top leg. It should be noted, however, that commensurate increases in drag force and power dissipation 
will also occur with the r-(f) design. Additionally, packaging constraints tend to limit the number of poles 
that may be wound in the r-(f> configuration to about one-half of the number of equivalent r-z poles. The 
physical and electrical constraints of the application of the axial Halbach magnetic bearing will usually 
determine the appropriate choice of stator orientation. 

These results were obtained using a small-scale test model. The factors limiting performance of this 
model actually improve dramatically when the physical scale is increased by an order of magnitude. First, 
the available space for packaging stator windings increases, permitting the use of lower gauge wire (less 
7?), a higher number of turns (more Z), or both. Increasing the radial span (W) of the stator winding gives 
a net gain in force production since Fi ift increases with W 2 , while the winding resistance only increases 
linearly with W. This is due to the fact that W figures linearly into both the stator current calculation and 
the Lorentz force calculation. For large machines on the order of 1 m diameter, the stator inductance to 
resistance ratio may be made favorable enough to make the use of an external series inductor 
unnecessary, even at speeds of less than 1000 rpm. This favorable balance results in less resistive heating 
and the larger size gives more thermal mass, as well as more surface area from which to radiate heat. 

Concluding Remarks 

The NASA Glenn Research Center has successfully designed, developed, analyzed, and tested a 
revolutionary Axial Halbach Magnetic Bearing. The goals of the project include improving aircraft 
efficiency, reliability, and safety. The objective of this work is to develop a viable non-contact magnetic 
thrust bearing utilizing Halbach arrays for all-electric flight, and many other applications. This concept 
will help to reduce harmful emissions, reduce the Nation’s dependence on fossil fuels and mitigate many 
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of the concerns and limitations encountered in conventional axial bearings such as bearing wear, leaks, 
seals, and friction loss. The Axial Halbach Magnetic Bearing is inherently stable and requires no active 
feedback control system or superconductivity as required in many magnetic bearing designs. The Axial 
Halbach Magnetic Bearing is useful for very high speed applications where conventional bearings cannot 
be used including turbines, instrumentation, and medical applications. 

Theoretical derivations have been developed successfully to predict the levitation forces generated by 
a circular Halbach array and coil assembly. Finite element analyses successfully validated the theoretical 
derivations. Empirical test results obtained from experimental hardware successfully validated the basic 
principles described, and the theoretical work that was performed. Of particular value, are the analytical 
tools and capability that were developed successfully under this project. Performance predictions can be 
made confidently for machines of various scale. 

The test results were obtained using a small scale test model. The factors limiting performance of this 
model improve significantly when the physical scale is increased. For large machines on the order of 1 m 
diameter, the stator inductance to resistance ratio may be made favorable enough to make the use of an 
external series inductor unnecessary, even at speeds of less than 1000 rpm. This favorable balance results 
in less resistive heating and the larger size gives more thermal mass, as well as more surface area from 
which to radiate heat. 

The report concludes that the implementation of Axial Halbach Magnetic Bearings can provide 
significant improvements in rotational system performance and reliability. In addition to aircraft engines, 
this technology has potential application in ultra-efficient motors, computer memory systems, 
instrumentation systems, medical systems, manufacturing equipment, and space power systems, such as 
generators and flywheels. 
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Appendix A 

Equipment Under Test Summary Data Sheet 


Axial Halbach Magnetic Bearing Test Model 


u 

Drive Motor 



1.1.1 Type 

Brushless DC Servo 


1.1.2 DC Voltage Range 

20 to 48 V 


1.1.3 Operating T emperature 

0 to 70 °C (32 to 158 °F) 


1.1.4 Storage T emperature 

-20 to 85 °C (-4 to 185 °F) 


1.1.5 Humidity 

0 to 80 percent 


1.1.6 Peak Torque 

0.63 N-m (89 oz-in.) 


1.1.7 Continuous T orque 

0.24 N-m (33 oz-in.) 


1.1.8 No Load Speed 

7400 rpm 


1.1.9 Nominal Power 

0.10 kW (0.14 HP) 


1.1.10 Servo Update 

4069 Hz 


1.1.11 Rotor Inertia 

1.34x10 5 kg-m 2 (0.0019 oz-i 


1.1.12 Width 

57.2 mm (2.25 in.) 


1.1.13 Length 

115 mm (4.5 in.) 


1.1.14 Weight 

0.95 kg (2.1 lb) 

1.2 

Force Measurement System 



1.1.2 Type 

Variable Capacitance 


1.1.3 Capacity 

410 g 


1.1.4 Readability 

0.001 g 


1.1.5 Linearity 

±0.02 g 

Rotor 

2.1 

Diameter 

11.94 cm (4.7 in.) 

2.2 

Weight 

660 g (1.46 lb) 

2.3 

Magnet Type 

Neodymium Iron Boron B55 

2.4 

Magnet Orientation 

Halbach Array 

2.5 

Magnet Shape 

Sectors 

Single Turn r-z Stator Pole Piece 


3.1 

Wire 

20 Square Copper 

3.2 

Inductance 

0.08 juH 

3.3 

Resistance 

7.55 m/2 
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4.0 

Six Turn r-z Stator Pole Piece 



4.1 

Wire 

20 Square Copper 


4.2 

Inductance 

1.90 juH 


4.3 

Resistance 

22.50 m/2 

5.0. 

Single 

Turn r-(f) Stator Pole Piece 



5.1 

Wire 

20 Square Copper 


5.2 

Inductance 

0.085 juH 


5.3 

Resistance 

10.00 m/2 

6.0 

Six Turn Stator Pole Piece 



5.1 

Wire 

20 Square Copper 


5.2 

Inductance 

2.34 nH 


5.3 

Resistance 

19.42 m/2 

7.0 

Inductor 



7.1 

Configuration 

Helix 


7.2 

Wire 

10 Square Copper 


7.3 

Number of Turns 

17 


7.4 

Core 

Ferrite Rod 


7.5 

Inductance 

9.1 juH 


7.6 

Resistance 

2.5 m/2 
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Appendix B 

System Performance Test Results 

A complete set of plots of the test results are included here. Table 1 identifies the tests that were 
conducted. 
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